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Abstract 

In this paper we improve the result about the polyconvexity of the energies from the family of isotropic 
volumetric-isochoric decoupled strain exponentiated Hencky energies defined in the hrst part of this series, 
i.e. 

f detF>0, 

W^AF) = { ^ 2fc 

[ +00 if det F < 0 , 

where F = V</5 is the gradient of deformation, U = %/ F'^F is the right stretch tensor and devnlogl/ is the 
deviatoric part of the strain tensor log U. The main result in this paper is that in plane elastostatics, i.e. 
for n — 2, the energies of this family are polyconvex for k > j, k > ^, extending a previous result which 
proves poly convexity for fc > |, fc > i. This leads immediately to an extension of the existence result. 


Key words: finite isotropic elasticity, Hencky strain, logarithmic strain, natural strain, polyconvexity, rank 
one convexity, volumetric-isochoric split, existence of minimizers, plane elastostatics, existence of minimizers. 


1 Introduction 

1.1 Exponentiated Hencky energy 

In a previous series of papers |29l 1301128] we have modified the Hencky energy and considered the family of 
energies 

„ ( I|dev2 logC/||^ fe [(logdetC/)]^ -r dpf F > f) 

lP,e(^) = T-(T4v)+lPT(deteKl)= (I.l) 

( +00 if detF<0. 

We have called this the exponential Hencky energy. Here, /i > 0 is the infinitesimal shear modulus, k = 
> 0 is the infinitesimal bulk modulus with A the hrst Lame constant, k, k are additional dimensionless 
parameters, F = is the gradient of deformation, U = V F'^F is the right stretch tensor and dev 2 logC/ = 
logC/ — itr(logl7) • I is the deviatoric part of the strain tensor logV. For X £ ||V|p = {X,X) is the 

Frobenius tensor norm, tr(X) = (X, 1) and 1 denotes the identity tensor on For further notations we 

refer to [H]. 
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Our renewed interest in the Hencky energy is motivated by a recent finding that the Hencky energy (not 
the logarithmic strain itself) exhibits a fundamental property. By purely differential geometric reasoning, in 
forthcoming papers EUllSlISI] (see also HI Ell it will be shown that 


( 1 . 2 ) 

where distgeod is the canonical left invariant geodesic distance on the Lie group GL“''(n) and distgeod,SL(n)j 
distgeod,R+.i denote the corresponding geodesic distances on the Lie groups SL(n) and M-j. • 1, respectively (see 
[13 [31]). On the other hand, some constitutive issues, e.g. the invertible true-stress-true-strain relation and the 
monotonicity of the Cauchy stress tensor as a function of log i?, where B = F is the left Cauchy-Green tensor, 
recommend the energies from our family of exponentiated energies as good candidates in the study of nonlinear 
deformations. Moreover, in the first part [H] it is shown that the proposed energies have some other very useful 
properties: analytical solutions are in agreement with Bell’s experimental data; planar pure Cauchy shear stress 
produces biaxial pure shear strain and the value 0.5 of Poisson’s ratio corresponds to exact incompressibility. 
It is found [Hj that the analytical expression of the pressure is in concordance with the classical Bridgman’s 
compression data for natural rubber. An immediate application to rubber-like materials is proposed in [24] . 
We have also shown that the energies from the family of exponentiated energies improve several features of the 
formulation with respect to mathematical issues regarding well-posedness. We have established that, in planar 
elasto-statics, the exponentiated energies satisfy the Legendre-Hadamard condition (rank-one convexity) 
[29] for k > k > i, while they are polyconvex [30] for k > ^, k > ^ and satisfy a coercivity estimate [30] for 
k > 0, k > 0. These results now allow us to show the existence of minimizers [30] for fc > i, k > i. 

1.2 Poly convexity 

The notion of polyconvexity has been introduced into the framework of elasticity by John Ball in his seminal 
paper [31 El [33]. Various nonlinear issues, results and extensive references are collected in Dacorogna [10] . In 
the two dimensional case, a free energy function W(F) is called polyconvex if and only if it is expressible in the 
form W{F) = P{F, det F), P : —>■ R, where F(-,-) is convex. Polyconvexity is the cornerstone notion for 

a proof of the existence of minimizers by the direct methods of the calculus of variations for energy functions 
satisfying no polynomial growth conditions, which is the case in nonlinear elasticity since one has the natural 
requirement W{F) —?► oo as det F —>■ 0. Polyconvexity is best understood for isotropic energy functions, but 
it is not restricted to isotropic response. The polyconvexity condition in the case of space dimension 2 was 
conclusively discussed by Rosakis [35] and Silhavy [42] 1331 Ell [33 [33 ES [44], while the case of arbitrary spatial 
dimension was studied by Mielke [21]. The n-dimensional case of the theorem established by Ball la page 
367] has been reconsidered by Dacorogna and Marcellini [T^, Dacorogna and Koshigoe [TT] and Dacorogna and 
Marechal m- It was a long standing open question how to extend the notion of polyconvexity in a meaningful 
way to anisotropic materials [3]. An answer has been provided in a series of papers [33[33[37|[I3[3a[a[3Q]. 
Rank-one convexity domains for the Hencky energy 

W^iF) = W„{U) = /r||dev„logC/f+ ^[tr(logt/)]2 

have been established in [53 E313 [IS]- Satisfaction of the Baker-Ericksen inequalities in terms of the logarithmic 
strain tensor is discussed in [3 [52], while necessary conditions for Legendre-Hadamard ellipticity are given in 

m- 

1.3 Motivation 

We have remarked that the rank-one convexity holds true for k > while the polyconvexity holds true for 
> i. Hence, the following question arose: is there a gap for ^ > /c > -j? In a previous work we have used 
the sufficiency condition for polyconvexity which has been discovered by Steigmann [491150] . Eventually, it is 
based on a polyconvexity criterion of Ball [3 , but it allows one to express polyconvexity directly in terms of the 
principal isotropic invariants of the right stretch tensor U (see also [131 [33 [13 [13 [13 13 E])- As it turns out, 


distgeod ((detF)^/” • l,SO(n)) = distg^e^i.R+.i ((detF)i/” -1,1) = jlogdetFp, 
‘iistgeod ( (detF)i/n ’^^(^)) = d'st2^od,SL(n) ( (detF) = l|dev„logl7f 
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in plane elastostatics, Steigmann’s criterion is already hidden in another sufficiency criterion for polyconvexity 
given earlier by Rosakis [36]. However, Steigmann’s criterion is clearly not necessary for polyconvexity. Hence, 
our previous results may be improved. 

In this paper, we use a direct approach based on the fact that the function Y : [l,oo) —>■ R given by 
Y{d) = 6 2 0 € [l,oo), is convex for fc > | and it is also increasing for 9 >1, while Z : GL+(2) —>• [l,oo) 

given by Z{F) = F G GL^(2), where Amax is the largest singular values (principal stretches) of F, is 

polyconvex. Therefore, we prove that the exponentiated Hencky energies are in fact polyconvex for k > j and 
the main existence results is also valid for these values of the fitting parameter k. 


2 Preliminaries. Auxiliary results 


2.1 Formulation of the static problem in the planar case 

The static problem in the planar case consists in finding the solution tp of the equilibrium equation 

0 = DivS'i(V<p) in (2.1) 

where the first Piola-Kirchhoff stress tensor corresponding to the energy W,,jj [F) is given by the constitutive 
equation 


Si (F) = [2/r e'”" log £^ll\ dev 2 \ogU + K u)f j 


F- 


e fi, 


( 2 . 2 ) 


with F = V(p, U = V F'^F. The above system of equations is supplemented, in the case of the mixed problem, 
by the boundary conditions 


(p{x) = 0i{x) on Td, Si{x).n ='si{x) on Tjv, (2-3) 

where T^i, Tjv are subsets of the boundary 9H, so that U T^v = dSl, r£inrjv = 0, nis the unit outward 
normal to the boundary and are prescribed fields. 


2.2 Auxiliary results 

In [29| and [30] the following results were established. 

Lemma 2.1. |29j Let fc G K and the matrix F G GL''’(2) with singular values Ai, A 2 . Then 


Wl^°{F) = e'=ll'iov 2 iogGf ^ =g(Ai,A 2 ), where g : ^ R, 5 (Ai,A 2 ) := 

Lemma 2.2. [29] Let m G N. Then the function 1 1 —>■ is convex if and only if k > — 

This lemma together with a results established in [B page 213] led to: 

Proposition 2.3. (Convexity of the volnmetric part) [2^ The function 

F ^ Wf°\F) := p ^ ql +(„) 


(log 


>1 > 


1 

GX+TT- 


2 


(2.4) 


is rank-one convex in F for k > ■ (More explicitly, for m = 2 this means k > ^, in case of m = 3 

rank-one convexity holds for fc > 

Theorem 2.4. (Coercivity) |30j Assume for the elastic moduli p > 0, k > 0 and k > 0, k > 0. Consider the 
energy I{(p) = W^^(V(p(x)) dx, where W^^{F) = W^^{U) = 2 log uW^ _|_ ^ gfc|tr(iogC/)p_ Then L{ip) is 

g-coercive for any I < q < 00 . 
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3 Improved polyconvexity result 




ssTT' for each 


By Lemma El the function W^^{F) = e''Ildev 2 log c/|| jg W^^°{F) = e 

F € GL+(2), where Ai > A 2 is an ordered pair of singular values of F. We view IV as the composition 
= Y o Z, where Y : [l,oo) —>• R is given by Y{9) = £ [l,oo), and Z : GL'''(2) —>• [l,oo) is given 

by Z{F) = , F £ GL“'"(2). For k> j the function Y{9) is convex and it is also increasing, see Lemma 

and Figure[T]. Hence, in order to prove the polyconvexity of fW™, in view of Lemmait is suffices to prove 



Figure 1: The graphic of the function Y : [l,oo) R, Y{d) = G [ 1 , 00 ), for 

different values of the fitting parameter k. 


that Z is polyconvex, since we have the following result: 

Lemma 3.1. If Y : [l,oo) —>■ R is increasing and convex and Z : GL''’(2) —>■ [l,oo) is polyconvex, then Y o Z 
is polyconvex. 

Proof. This is proved in exactly the same way as |341 Theorem 5.1]; only the elements x and y of the proof have 
to be interpreted as pairs {F,S), where F £ GL“''(2) and A > 0. The generalization to arbitrary dimension, if 
needed, is straightforward. □ 

Lemma 3.2. The function Z is polyconvex on GL^(2). 

Proof. An elementary evaluation of the Hessian matrix shows that for p > 0, g > 0 the function /(a, b) = ^ is 
convex on [0, 00 ) x [0, 00 ) if and only if g < p — 1. We take p = 2 and q = 1 and we observe that 

Z(F) =/(Ai,detV), (3.1) 

for every F £ GL“''(2). Let fi{F,S) = f{\i,S), for each F £ GL^(2) and ^ > 0. Using the well-known fact 

that F I—>■ Ai(F) (the largest singular value) is convex (see e.g. [31 Lemma 5.3] and also Annendix lA.ll) . we now 

prove that fi is convex. Indeed, since / is increasing in the first argument, we have for any t, 0 < t < 1, and 
any Fa £ GL'''(2), a = 1,2, that 

/i((l -t)Fi+tF 2 ,{l-t)Si+t 62 ) = /(Ai((l -t)Fi+t F 2 ), {l-t)5i + t 82 ) (3.2) 

< (1 ~ ^) /(^i(Vi), (5i) -h t f{Xi{F 2 ), 82 ) = (1 — t) fi(Fi, 8 i) + t fi{F 2 , 82 ), 

where we used the convexity of / asserted above. Thus Z{F) = fi{F,det F) is polyconvex. □ 

Therefore, we conclude: 

Proposition 3.3. If k> j, then the function F 1 —lFg°(U) = IIiog£||^ jg polyconvex on GL“'"(2). 

In view of the above proposition and using Corollary 12.31 we conclude that: 

Theorem 3.4. The functions W],jj : R^^^ —>■ R from the family of exponentiated Hencky type energies are 
polyconvex for p>0,K>0, fc> — and k > —. 
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4 Existence result 


In plane elastostatics, having proved the coercivity and the polyconvexity of the energy W^^{U) for fc > | and 
fc > |, it is a standard matter to improve the existence results established in [30]. 

Theorem 4.1. (Existence of minimizers) Let the reference configuration C be a hounded smooth domain 
and let Tu be a non-empty and relatively open part of the boundary dLl. Assume that I{ip) = Jq Lp{x))dx 

where W^^{F) = W^^{U) = f ^ gfc |tr(iog£/)|^_ g q>l be given with I{ipo) < oo 

and p, > 0, K > 0, k > j and A: > i. Then the problem 

min|/(<p)= J W^^{Vip{x))dx, (p{x) = (Pq{x) for x G Td C dLt, G (4.3) 

admits at least one solution. Moreover, ip G W^'‘^(Ll), q > 1. Since we do not know whether det Vt/J > c > 0 
and since the energy does not satisfy a polynomial growth condition, it is not clear whether the Euler-Lagrange 
equations are satisfied in a weak sense. 
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Appendix 

A.l The convexity of the largest singular value of F 

In this section we outline the proof of the convexity of the function F hA Ainax(^) (the largest singular value of F). First, we 
present three lemmas given in [S] page 364]. 

Lemma A.l. (von Neumann El; see also Mirsky [221123j ^ Let A, B £ have singular values ai > cx 2 > ••• > c^n > 0 and 

01 ^ 1^2 > ••• > > 0. Then |tr(AB)| < (o, 0). 

Lemma A. 2. (von Neumann El; see also Mirsky [221123] ) Under the hunotheses of Lemma \ A. 1\ max \{AQ, B'^)\ = {ol, ^). 

Q,J?eO(n) 

Proof. There exist orthogonal matrices Qi<,Q 2 i Rit R 2 such that A = Qi diag(Q) i?i, B = Q 2 diag(/3) i? 2 - Choose Q = 

R = R 2 Q 1 . Then tT{AQ B R) = {a,l3). But for any orthogonal Q,R the matrices AQ and BR have singular values a, 0 
respectively. Hence \{AQ,R^ B'^)\ = tr^AQ B R) < {ol,0) by Lemma |A.II □ 

Using the above two lemmas, Ball (3] page 364] proved: 

Proposition A.3. (Ball [3] page 364] Let r\ > r 2 '> ■■■ > Tn > 0. Then (r. A) is a convex function of F, where Ai > A 2 > ... > 
An > 0 are the singular values of F. 

Proof. In Lemma |A.21 put A = F and B = diag(ri, r 2 ,..., rn). Then (r, A) = maxg tr{F Q B R). Since each tr(F Q B R) is 
a convex function of F, we obtain that (r. A) is also a convex function of F. □ 

Moreover, it holds true that: 

Corollary A. 4. The function F 1 —>• Amax = Ai = maxi=i^...^n Ai(F) is a convex function of F, where Ai > A 2 > ... > An > 0 are 
the singular values of F. 

Proof. By letting r = (1, 0, 0,..., 0) in Lemma lA.3l it follows that for Ai(F) is a convex function of F. □ 
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